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SYMMETRIC PRIMAL-DUAL SIMPLEX ALGORITHM 
FOR LINEAR PROGRAMMING 


Keshava Prasad Halemane 
Dept. of Mathematical and Computational Sciences, 
Karnataka Regional Engineering College, Surathkal, 
P.O. Srinivasanagar - 574 157 (D.K) 
Karnataka State, India 


ABSTRACT: 

A novel iterative algorithm to solve linear programming (LP) problems based upon 
a newly identified set of four (and/or possibly six) distinct types of simplex pivoting steps 
defined over a symmetric primal-dual pair of LP, is presented. The two (or possibly four) 
types of classical (standard) simplex pivots are the Primal Standard Pivot with positive 
(and/or zero) indicator, and the Dual Standard Pivot with negative (and/or zero) indicator. 
The newly identified pivot types are: the Primal Tricky Pivot with positive indicator and 
the Dual Tricky Pivot with negative indicator. 


The Symmetric Primal-Dual Simplex Algorithm is an iterative solution strategy, wherein 
each iteration consists of the selection and application of exactly one pivoting step, chosen 
in a well defined systematic preference order (for our convenience, at least), from among the 
possibly four (or even six) distinct types of pivots, that may be valid depending on. the 
tableau data entries. , 

For expressional convenience, this algorithm is discussed using the Tucker’s Compact 
Simplex Tableau for Linear Programming problems expressed in the standard (canonical) 
form representing the Symmetric Primal-Dual Pair. An analysis of the evolution of the tableau 
data entries pattern observed as iterations proceed, is also presented. A classification of the 
LP problem into one of the six types is identified based on the data entries pattern observed 
in the final tableau. 


1. INTRODUCTION: 
Linear Programming (LP) problem represents one of the most widely. used class of 


numerical/quantitative computational models, for which any possible improved solution. 


technique would certainly be highly desirable. Of course, there has been several alternative 
solution strategies suggested including the classical simplex method of Dantzig [1] and several 
variations thereof, followed by recent polynomial time algorithms, namely the Ellipsoid 
Method of Khachiyan [2, 3] and the Karmarkar Algorithm [4] both classified now as belonging 
to Interior Point Algorithms. 


Gn this paper, a novel generalization of the classical Simplex Method of Dantzig is 
presented. The fundamental concept in Dantzig’s Simplex Method consists of moving from 


one basis (tableau) to a neighbouring basis (tableau), by a single exchange between a non-basic . 


(entering) variable and a basic (leaving) variable, appropriately selected. We associate the 
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term "Simplex Algorithm" for our method because this fundamental feature is being preserved. 
The Proposed algorithm is an enhancement over the Dantzig’s Simplex Method, in terms of 
providing a wider scope for the selection of the pivots, as can be seen in the ensuing 
discussions. The development of the algorithm is mainly supported by the identification of 
a set of possibly four (or even six) distinct types of simplex pivot selections, maintaining a 
sense of symmetry between the primal and the dual problems. Also, for expressional efficiency 
or at least convenience, the solution Strategy and the analysis thereof are all based on the 
Tucker's Compact Simplex Tableau for LP expressed in the standard (canonical) form 
representing the symmetric primal-dual pair. Hence we justify the chosen name for the present 
algorithm: "Symmetric Primal-Dual Simplex Algorithm for Linear Programming”. 

2. PRELIMINARIES: 


We will go through some well known preliminaries for the sake of establishing the 
notational conventions used in this paper. The symmetric primal-dual pair of LP in the 
standard form is as follows. 


Primal Problem: 


Maximize cx =f 
s.t. : A.x <b (1) 
0 
Dual Problem: 
Minimize v.b=g 
S.t. vA2c (2) 
v20 


The matrix-dimensional descriptions for each of the problem parameters/data in (1) and 
(2) above are as follows: 


x Primal decision Variables nx! vector 
c Primal objective function coeffs 1xn vector 
f Primal objective function value 1x1 scalar 
A Primal constraint coeffs matrix mxn matrix 
b Primal constraint upper. bound mx vector 
v Dual decision variables Ixm vector 
g 


Dual objective function value 1x1 scalar 


We introduce the mx! vector y of slack variables to (1) and the 1xn vector u of surplus 
variables to (2) to write the symmetric primal-dual pair in the canonical form as follows 


Primal Problem: 
Maximize cx+Oy=f 
st. Ax +Imy =b (3) 
x,y2 0 
Dual Problem: 
Minimize v,b +u.0 =B 
s.t. v.A-u.ln = ¢ 
v,u20 


(4) 
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This primal-dual pair is represented in Tucker's Compact Simplex Tableau as: 


=e (5) 


-l 


I 
uj g 


For the LP problem pair (1) & (2) or equivalently (3) & (4) the tableau (5) represents 


the initial tableau indicating the initial basic solution (IBS) wherein y; are the primal basic 


variables in the initial basis associated (one to one permanent association) with v; the dual 
non-basic variables. and x; are the primal non-basic variables associated (one to one permanent 
association) with uj the dual basic variables, in the initial basis. Note that xj (and a -!) are 
column-labels and v,; (and a -i) are row labels in the tableau (5), and the way to interpret 


(read) the tableau is as follows 
Primal Problem: 
¥ aj -X 7 b,; = — y; i € R (row index) (6) 


jec 


¥ c -X — 0 = f (to be maximized ) 
jec 
Dual Problem: 
> yay G = Ye j € C (Column index) (7) 


eR 


. vy, -b;- 9 = g (to be minimized) 
sER 
wherein the variables x;. yj. Vj, Uj are all considered to be non-negative 
3. ALGEBRA (ARITHMETIC) OF THE PIVOTING PROCESS: 
It is an interesting 
selected, the actual pivoting process (the algebra an 
same irrespective of the p 


Pp 
pivoling process. 


For the sake of generality, let us imagine that we are somewhere in the middle of 


solving a LP problem, and have the system mode] represented by a tableau as follows: 
16 


but simple enough exercise to observe that once a pivot element is 
d hence the arithmetic operations) is the 
ivot selection; for example whether it is a primal pivot or a dual 
ivot. Hence it suffices to present here a single (common) set of operations representing that 
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, (8) 


By the nature of the sequence of elementary row (column) operations being performed 
during any pivoting process, the system model represented by (8) is equivalent to that 
represented by the initial tableau (5) which corresponds to the primal-dual pair (6) & (7). 
The transformed version of the primal-dual pair directly expressed by (8) is as follows 


Primal: 
B : 
z; = 8B; - > age Zh i € R(row index) 
jec 
f= - 6+ * 4-3 (to be maximized) (9) 
: jeC 
Dual: 
wP =-¥+ ¥ wh. a; j € C (column index) 
ieR 
g =-d5+ by wh . Bi (to be minimized) (10) 
ieR : 


The effect of a pivoting operation on (9) & (10) performed with a chosen pivot element 


Clyy is exactly to affect an exchange between the variables indicated by 1 and J in (9) and 
2B in (9), and w)* is entered 


(10). That is, Zz is entered into primal basis in exchange for 2; 
in (LO), Suppose we have chosen the pivot element 


into dual basis in exchange for wf 
ection scheme, and we would like to derive the resulting 


(yy Using some appropriate pivot sele 
tableau. Let the resulting tableau be indicated as: 
19 
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~ (By’ 


(11) 


(f)’ 


Then the process of deriving (11) from (8) constitute the algebra (arithmetic) of the 
pivoting process and is detailed below (without any derivational justifications): 


(yy) —(1/ayy): (a) <—  (a))/ OF; : (B)’ <— (By/a43): 
(Qi3)’ — - (G53! Oy); (yy) — — (Yy/Oy3); 

(BY — B - (Byoy) oy 
(5)’ — 8 - (Boys) vy 5 


followed by an exchange of the label associated with the row I and column J; that is 
effectively: 


(aj) — Oj — (Gy / ay) ay: 


(yi) — Yj - (yj / Oy) yy 


BY < 2: (PY zy; 

Gye a (zBY zB; 
and 

(apy owes (why we 

(wey — we (why wh; 


fori ER \ {I} andj EC \ {J} 


This completes the pivoting process. 


4. SIMPLEX PIVOT SELECTION SCHEMES: 
There are four (two pairs) fundamental types of pivot selection schemes namely Primal 
Standard Pivot (PSP), Dual Standard Pivot (DSP), Primal Tricky Pivot (PTP) and Dual Tricky 


Pivot (DTP). Further a primal/dual standard pivot can again be classified into one with’ 


positive/negative indicator (that is PSPPI and DSPNI) and one with zero indicator (that is 
PSPZI and DSPZI), thus resulting in a set of possible six distinct types of pivot selection 
schemes, that are available for simplex pivoting process in solving linear programming 
problems. The algebra of the pivot selection schemes are given below, along with a schematic 
representation of the tableau data entries pattern that leads to such pivot selection. 

20 a 
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4.1 Dual Standard Pivot with Negative Indicator, DSPNI: 


Je arg min ((y; / oy) | Bp <0; ¥, < 0; a; < 0} 
Ea re 


Je{jec ly,>0); 
I © arg min ((B, / a3) |B =0; y; > 0; ajy > 0} 
icgR 


Je(jec ly>0}; 


max 
I — arg se ((B, / a) | B <0; yy > 0; a; < 0} 
ie€R 


le{ieR |B <0}; 


J <— arg max {(7j / ty;) | By <'03 ¥> 0; ay; > 0} 
jec 


4.5 Dual Standard Pivot with Zero Indicator, DSPZ1: 
0 le le({ieR |B =0}; 
if J © arg min {(y) / ay) | B, = 0: Yj S 0; ay < 0} 
jec 7 


L+loele] 
L 
4.6 Primal Standard Pivot with Zero Indicator, PSPZI: 


Je{jec ly, =}: 


1 & arg min {(B / ay) | B 20: yy = 0; ay > 0} 
igR 


It can be easily observed that a primal (dual) tricky pivot with zero indicator is essentially 
the same as dual (primal) standard pivot with negative (positive) indicator, wherein the 
resultant minimum-ratio comes out to be zero, 


21 
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5. EFFECT OF PIVOTING OPERATION: 


It is useful at this point to make a few observations regarding the effect of pivoting 
operation. in each of the above pivot selection schemes. 


DSPNI brings about an immediate improvement in the primal feasibility (w.r.t. the pivot 
row, at least) without deterioration of dual feasibility. The extent of this improvement in 
primal feasibility can be measured by the corresponding resultant improvement (decrease) 
observed in the value of the dual objective function, given by | By; /a1) | - 


PSPPI brings about an immediate improvement in the dual feasibility (w.r.t. the pivot 
column, at least) without deterioration of primal feasibility. The extent of this improvement 
in dual feasibility can be measured by the corresponding resultant improvement (increase) 
observed in the value of the primal objective function, again given by | Byyy /Qyy | - 


PTPPI brings about an immediate improvement in the primal feasibility (w.r.t. the pivot 
row. at least) without any concern to the dual feasibility. The extent of this improvement in 
primal feasibility can be measured by the corresponding resultant improvement (deerease) in 
the value of the des, objective function, again given by | By; /oy | - 

L 


DTPNI brings about an immediate improvement in the dual feasibility (w.r.t. the pivot 
column, at least) without any concer to the primal feasibility. The extent of this improvement 
in dual feasibility can be,measured by the corresponding resultant improvement eee 
the value of the Frama objective function, again given by | By; /o, | - 


The other two types of pivot selections namely PSPZI and DSPZI are utilized to affect 
a change in the basis, in situations with multiple-solutions and/or degeneracy, and do not in 
any way affect the primal or dual feasibilities. However they certainly play an important rolé 
in the overall] solution strategy. 


From the above, one can observe that as a general, uniform and common basis for the 
local effectiveness measure of a particular pivoting operation (applicable for any and every 
iteration. for both primal] and dual) we can use the absolute value of the change in the value 
ef the objective function, given by lem(I,J) = abs(Byyy /O,)). It is not specifically suggested 
here, but one can ept to choose a pivot to maximize this local effectiveness measure in every 
ileration either among the possible pivots of a particular type or even among of all the 
possible pivots of all the possible types. Even if done so, it cannot be guaranteed (as per 
worst case analysis) 10 minimize the overal] number of iterations required for reaching an 
optimum solution. lt requires further research work to thoroughly understand, analyze and 
incorporate the concept of “local effectiveness measure (lem) of a single simplex pivoting 
operation” to the fullest extent, in developing an efficient solution strategy. For now, let us 
come to the main algorithm itself. 


6. SYMMETRIC PRIMAL-DUAL SIMPLEX ALGORITHM: 


The six distinct types of pivot selection schemes are considered in the following order 
of preference: 


{{DSPNI, PSPPI}, {PTPPI, DTPNJ), {DSPZI, PSPZI}} 
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At every iteration, an attempt is made to select a pivot element, by checking the possible 
pivot selections belonging to one of the above six types of pivot selection schemes in the 


preference order pre-specified. 
That is, a pivoting operation is executed with a pivot element of the type- 
i. DSPNI if possible, 
ii. _PSPPI only if none of DSPNI is possible, 
iii. PTPPI only if none of DSPNI nor PSPPI is possible, 
iv. DTPNI only if none of DSPNI, PSPPI, PTPPI is possible, 
and when none of the above four is possible, go for 
v. DSPZI and/or 
vi. PSPZI in order to explore the possible alternative basic solutions. 


It is to be noted that depending upon the actual data entries in the tableau, a pivot 
selection of specific type which was not possible in an earlier iteration, can become possible 
in a later iteration, sometimes even in the very next following iteration. That is why it is a 
crucial part of the algorithm to check in each and every iteration, for each of the six types 
of possible pivot selection schemes in the pre-specified order. 


7. EVOLUTION OF TABLEAU DATA ENTRIES PATTERN: 


As a means to show the convergence of the algorithm to one of the possible six 
termination patterns (discussed later), an analysis of the evolution of the tableau data entries 
pattern is given here. Figure-1 gives the six different patterns applicable when each of the 


Figur 1. Evolution of tableau data entries pattern. 
- negative O Zero + positive 

0 non - positive * any value @& non-negative 
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six different pivot selection schemes may be chosen as a possibility in some iteration during 
the process of solving a LP problem. This figure needs to be read keeping in mind the algebra 
of the pivoting process and the details regarding the nature of the pivot selection itself. To 
facilitate the analysis of the possible termination patterns, it is useful to look at a systematic 
classification of all the 27 possible combinations of tableau data entries pattern aS shown in 
Figure-2. Here. the possible situations that would lead to a successful pivoting are indicated- 
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Figure 2. Twenty seven distinct combinations of tableau data entries pattern, 
~ negative O Zero + positive 
* any value 
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they add up to ten, including the four situations corresponding to multiplicity/degeneracy. 
Others wherefor no pivot selection is possible, are marked according to the nature of the 
termination, by indicating F for basic-feasible-finite, ® for infeasible and - for non- 
basic-feasible-infinite, corresponding to both primal and dual variables. Note that this 
classification is as per the quite well established approach, although a refinement is indicated 
to distinguish the case wherein primal (dual) has an infinite non-basic optimum with finite 
value for the objective function while the dual (primal) has a feasible-finite optimum. This 
distinction from the classical approach arises because we give primary emphasis on the 
classification based on the nature of the decision variables at termination, and give secondary 
emphasis on the finiteness (or otherwise) of the objective function value. 


Thus the seventeen (out of the total of twenty seven) situations that lead to non pivoting 
termination can be classified into six distinct classes as follows 


i) Six situations with P = F, D= F 
ii) One situation with P =F, D= oe 
iii) One situation with P = ‘co; D=F 
iv) Four situations with P = ~, D=®90 
v) Four'situations with P=9®, D=o0 


vi) One situation with - P= 60, D=®@® 
Four among the remaining ten (out of the total of twenty seven) situations that indicate 
pivoting possibilities can be included in case (i) above and treated as possible termination 
situations because further pivoting would be useful only to the extent of discovering multiple 
optima for the primal/dual. The other six correspond to the six types of pivot selection 
schemes indicated. Figure-3 summarises these observations. 


From the above discussion and the observations presented in Figures 1,2 and 3 it is 
clear that the algorithm converges to one of the six distinct classes of termination situations 
and the given LP-Problem can be classified accordingly. 


8. SOME GENERAL COMMENTS: 


It is worth noting here that using the above proposed algorithm one can solve any LP 
problem, by first converting it into the standard/canonical form and then applying the proposed 
algorithm. In performing such transformation, it is possible to enhance the overall efficiency 
by the following approach: 


1) Free variables can. be replaced by non-negative variales not double in number, but 
only just one extra in number 
ii) Equations can be replaced by inequalities not double in number, but only just one 
extra in number. 
iil) No need for use of artificial variables. 
iv) The initial basic solution need not necessarily be feasible. 
Degeneracy and possible cycling can be handled in the usual method, by lexicographic 
ordering of the variables (equivalent and in fact supported by the concept of polynomial 
perturbations technique). 
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P=oa 


Figur 3. Six categories for tableau data entries Pattern at termination. 
— negative O Zero + positive 
6 non-positive * any value © non-negative ¢ un-analyzed 


Also it is possible to have alternative ordering of the pivot selection schemes incorporated 
in the algorithm within each of the three subsets indicated in section 6 above, although the 
relative ordering of these three subsets themselves is not to be altered. 


Although rigorous mathematical analysis of the computational complexity of this algorithm 
is yet to be studied in detail, it seems very likely to perform at least as good as (if not better 
than) the existing solution techniques for linear programming. 
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GEOMETRIC GRAPHS 


S. M. HEGDE 
Department of Mathematical and computational Sciences 
Karnataka Regional Engineering College 
P.O. Srinivasanagar - 574 157. Karnataka 


ABSTRACT 
A (p,q)-graph G is said to be (a, r) -geometric if its vertices can be assigned distinct 
positive integers so that the values of the edges, obtained as the products of the numbers 
assigned to their end vertices, can be arranged in the geometric progression a, ar, 
ar’, 5 ar™”. In this paper we prove that some well known classes of graphs are geometric 
for certain values of a and r and also initiate a study on the structure of finite (a, r) - 
geometric graphs. 


1. INTRODUCTION 


For all standard terminology and notation in graph theory we refer Harary [3]. 


Several practical problems in real life situations have motivated the study of labelings 
of a graph G = (V, E), which are required to obey certain variety of conditions depending 
upon the structure of G such as adjacency. There is an enormous amount of literature built 
up on several kinds of labelings of graphs over the past three decades or so, and it would 
be too unwieldy to list them here. For a cursory look at the topic we refer the reader to 
Bloom [2]. \ 


Labeling of G is merely an assignment of certain distinct nonnegative integers to the 
vertices of G. : 


In this paper we are interested in the study of vertex functions f: V(G) > A, A CN. for 
which the induced edge function f*: E(G) — N is defined as f* (uv) = f(u)f(v), Vuv € E(G). 
Such vertex functions are said to be multiplicative and henceforth this particular map f of 
f will be denoted f*. 


THEOREM 1 (Acharya and Hegde [1]): For any graph G and for any multiplicative 
vertex function f : V (G) > N 
(1) Mm fe) = A fw 
ee E(G) ue V(G) 
COROLLARY 1.1: If G is an r-regular graph then for any multiplicative vertex 
function f of G, 


nm @(e) = Tl f(u)'. 
e€ E(G) ue V(G) 


2. MULTIPLICATIVE NUMBERINGS 
A multiplicative numbering of a graph G is an injective multiplicative vertex function 


{ such that the induced edge function f* is also injective. 
af 
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